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1. Introduction 



Riemann zeta function for Re(s) > 1 can be defined as follows: 

OO 1 

c(*) = E^ (i-i) 

n=l U 

It also can be represented as a product with respect to all prime numbers p 

as) = u^-p- s r 1 (i-2) 

p 

It can be analytically continued in the whole complex plane of s. It has only one pole, at 
s = 1. Riemann zeta function satisfies a functional equation 

C(s) = 2 V 1 sin(s7r/2)r(l - s)((l - s) (1.3) 

It has 'trivial' zeros at s = — 2n ( n > 1 is an integer). The famous Riemann hypothesis 
[15] states that nontrivial zeros belong to the straight line Re(s) = 1/2. Recently Montgomery 
and Odlyzko conjectured that for large values of imaginary part of s the distribution of zeros can 
be described by GUE of random matrices, see [25] and [24]. Forrester and Odlyzko related the 
problem of distribution of zeros to Painleve differential equation and integrable integral operators 
[21]. Riemann zeta function is useful for study of distribution of prime numbers on the real axis 
[14]. The values of Riemann zeta function at special points were studied in [17] , [18] . At even 
values of its argument zeta function can be expressed in terms of powers of n and Bernoulli's 
numbers 

C(2n) = (-l) n+l 2 2n -\ 2n B 2n /(2n)\ (1.4) 

The values of Riemann zeta function at odd arguments provide infinitely many different irrational 
numbers [16] . Riemann zeta function plays an important role, not only in pure mathematics but 
also theoretical physics. Some Feynman diagrams in quantum field theory can be expressed in 
terms of C( n )j see ? f° r example, [1] . In statistical mechanics Riemann zeta function was used 
for the description of chaotic systems [19]. One can find more information and citation on the 
following web-cite http://www.maths.ex.ac.uk/ mwatkins/ . 

We argue that ((n) is also important for exactly solvable models. The most famous integrable 
models is the Heisenberg XXX spin chain. This model was first suggested by Heisenberg [3] in 
1928 and solved by Bethe [4] in 1931 . Since that time it found multiple applications in solid state 
physics and statistical mechanics. 

The Hamiltonian of the XXX spin chain can be written like this 

H = J2 KX + i + ^f+i + - 1 ) (1-5) 

i=i 

Here N is the length of the lattice and of, of, of are Pauli matrices. We consider thermodynamics 
limit , when N goes to infinity. The sign in front of the Hamiltonian indicates that we are 
considering the anti-ferromagnetic case. We consider periodic boundary conditions. Notice that 
this Hamiltonian annihilates the ferromagnetic state [ all spins up]. 
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The construction of the anti-ferromagnetic ground state wave function | AFM > can be credited 
to Hulthen [5]. An important correlation function was defined in [8]. It was called the emptiness 
formation probability 

n 

P(n) =< AFM\ [] Pj\AFM > 
i=i 

where Pj = (1 + <t|)/2 is a projector on the state with spin up in jth lattice site. Averaging is over 
the anti-ferromagnetic ground state. It describes the probability of formation of a ferromagnetic 
string of the length n in the anti- ferromagnetic background \AFM > . In this paper we shall first 
study short strings ( n is small), in the end we shall discuss long distance asymptotic ( at finite 
temperature). The four first values of the emptiness-formation probability look as follows: 

P(l) = \ = 0.5, (1.6) 

P(2) = -(1 - ln2) = 0.102284273, (1.7) 

P(3) = \ - m2 + -C(3) = 0.007624158, (1.8) 
4 8 

P(4) = - - 21n2 + — C(3) - — C(3) hi2 - — ( 2 (3) 
55 85 

- 2|C(5) + 2|C(5) ln2 = 0.000206270 (1.9) 

Let us comment. The value of P(l) is evident from the symmetry, P(2) can be extracted from 
the explicit expression of the ground state energy [5]. P(3) can be extracted from the results of 
M.Takahashi [9] on the calculation of the nearest neighbor correlation. It was confirmed in paper 
[10]. One should also mention independent calculation of P(3) in [11]. One can express P(3) in 
terms of next to the nearest neighbor correlation 

< S*S* +2 >= 2P(3) - 2P(2) + lp(l) (1.10) 

The calculation of P(3) and P(4) is discussed in this paper. 
The expression above for P(4) is our main result here. 

We briefly announced our results in [2], here we provide the detailed derivation. The plan of 
the paper is as follows. In the next section we discuss a general procedure of the calculation of 
P{n). We also show how this scheme works for P(2). In the Appendices A and B we describe in 
detail the calculation of P(3) and P(4) respectively by means of the technique elaborated in the 
Section 2. The main results are summarized in the conclusion. 

2. General discussion of the calculation of P(n) 

There are several different approaches to investigate P(n): 

• Representation of correlation functions as determinants of Fredholm integral 
operators .This approach is based on following steps: 
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i. Quantum correlation function should be represented as a determinant of a Fredholm 
integral operators of a special type. We call these operators integrable integral operators. 

ii. The determinant can be described by completely integrable equation of Painleve type. 

iii. Asymptotic of correlation function [and the determinant] can be described by Riemann- 
Hilbert problem. 

This approach was discovered in [23] , it is described in detail in the book [6] . 

It is interesting to mention that this approach was successfully applied also to matrix models 
[22] 

• Vertex operator approach was developed in Kyoto by Foda, Jimbo, Miki, Miwa and 
Nakayashiki. This approach is based on study of representations of infinite dimensional 
quantum group U q SL(2), see [7] . 

We shall use the integral representation obtained in [8] 

p( n ) = f dXl f dX * f dXn f\(l I 1 )n-a( nX * )n TT SUlll 7r(Aj — Afc) 

1 ; Jc 2mA! Jc 2m\ 2 ' " Jc 2mX n fi/ Xj l sinh7rA/ ^ n n(Xj - X k - i) ' 1 ' ' 

The contour C in each integral goes parallel to the real axis with the imaginary part between 
and —%. In the frame of algebraic Bethe Ansatz this formula was derived in [20]. Recently such 
formula was generalized in paper [12] to the case, where averaging is done over arbitrary Bethe 
state [ with no strings ] instead of anti-ferromagnetic state. 

Let us describe in general a strategy that may be used for the calculation of P(n). The integral 
formula (2.1) can be easily represented as follows: 

n r d\ 

P ( n ) = Iij c ^ U n (\ u A n ) T n (X u . . . , An) (2.2) 

r . ( A,,....A.) = JB^q^TL ( , 4) 

First of all, let us note that in principle the contour C can be chosen between and —i 
arbitrary. Let us denote C a the contour that goes from ia — oo to ia + oo. In what follows it will 
be convenient to choose a = — 1/2 i.e. to integrate over the contour C_i/2- 

As appeared we can make a lot of simplifications without taking integrals but using some 
simple observations and properties of the function in the r.h.s. of (2.2) which has to be integrated. 

Let us define a "weak" equality ~. Namely, let us say that two functions F n (Ai, . . . , A n ) and 
G n (Xi, . . . , A n ) are "weakly" equivalent 

F n (A 1; ...,A n ) ~ GV t (Ai, . . • , A n ) (2.5) 



where 



and 
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if 

dx 



II / ^ U n (X u ...,\ n ) (F n (Ai, . . . , A n ) - G n (\ u . . . , A n )) = 0. (2.6) 

j=l ^ C-i/2 

Let us also introduce a " canonical" form of the function by the following formula 



^(Ai, ...,X n ) = ± if fl (2.7) 

where Pj n ^ are some polynomials of the form 

Pj ^ = -P,- ^(Ai, A3, . . . , A 2 j_i|A 2 j+i, A 2j ,+2, • • • , A n ) = 

U U, i 3 ,...,i 2j --i A l A 3 ••• A 2j-1 A 2j+1 A 2j+2 ••■ K 

< ii , 13, ■ 

< i 2j + i < i 2i 

where 



u ^ l\ , 23 , . . . , l2j — 1 ^ n — ^ 
< i2j + l < 4 2j+2 < ■ ■ ■ < i n < n — 1 



^2j+l,l2j+2,---,ln -/3 / £ r*2j + l,«2j+2, -,'n 

with j3 = or 1 in accordance with the equality 

P + H+13 + ■ ■ ■ + i 2 j-i + «2j+i + i2j+2 + . . . + i n = j + n mod 2 

and some rational numbers Cj 23+1 '? 23+2 'i''' 1 " . 

This form has some arbitrariness because if we substitute Xj = Xj — i/2 where all Xj are real 
then it is easy to see that the function U n (xi, . . . ,x n ) = U n (x 1 — i/2, . . . ,x n — i/2) transforms 
when {xi, . . . ,x n } — * {—x±, . . . , — x n } as follows 

~ n{n — 1) ~ 

U n (-x 1 ,...,-x n ) = (-l)^ - U n (x u . . . ,x n ). (2.9) 
Therefore any function F n (x±, . . . , x n ) that satisfies 

F n (- Xl , -x n ) = (-1)t+i p n ( Xl , ...,x n ) (2.10) 
being integrated makes zero contribution 



II / t^^i,---,^) F n (x 1 ,...,x n ) = 0. 

j = l J —00 Z7Tt 



It means that in order to get a nonzero result one should have the function F n {x\, . . . , x n ) of the 
same parity as of the function U n (xi, . . . ,x n ). Then if we re-expand the form (2.7) in terms of 
variables Xj instead of Xj we can fix the arbitrariness by imposing some additional constraints, 
namely, 

Pj ^3; • • • ? x 2j-l\%2j+l, x 2j+2, ■ ■ ■ , %n) = 

= P < f l \x 1 - i/2,x 3 - i/2, . . .,x 2 j-i - i/2\x 2 j+i - i/2,x 2j+2 - i/2, ...,x n - i/2) = 
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»l i i 3 , »2j-l X l X 3 • • • X 2j-1 x 2j+l x 2j+2 ■ ■ ■ x n ■ 

< »l, <3, ■ ■ ■ , <2j-l < " — 2 
< *2j + l < ^2j + 2 < ■ ■ ■ < in < ^ — 1 
»1 + <3 + • • • + 4 2j-2 + 123 + 1 + *2j + 2 + • • • + i n = 3 + 11 mod 2 

(2.H) 

In comparison with the coefficients C'^ +1 ^^'- 2 '^_ n 1 which can be pure imaginary all the coefficients 

C^ 3+1, ^ 3+2 are real and rational numbers. 
So we can expect that the function 

Tftxi, . . • , x n ) = T„ c (x! - i/2, ...,x n - i/2) (2.12) 

should satisfy the following property 

TZ(- Xl , -x n ) = (-I)t fZ( Xl , ...,x n ) (2.13) 

Below the property (2.13) will be implied when we will speak about the "canonical" form (2.7-2.8). 
Besides, one can note that the function T°(x\, . . . , x n ) should be real for real variables Xj. 

Our hypothesis is that for any n one can reduce the function T n defined by (2.4) to the canonical 
form i.e. there exist polynomials Pj in (2.7) such that 

T n (\i : ■ ■ ■ , A n ) ~ T£(\ u . . . , X n ). (2.14) 

Unfortunately, for the moment we do not have a proof of this statement for any n but we will 
demonstrate below how it works for n — 2, 3, 4. 

In fact, the problem of the calculation of P(n) given by the integral (2.2) can be reduced to 
the two steps. The first step corresponds to the obtaining of the "canonical" form for T n . The 
second step is the calculation of the integral by means of this "canonical" form. 

To do this one needs the following simple facts: 

I. Since the function U n (Xi, . . . , A n ) is antisymmetric with respect to transposition of any 
pair of integration variables, say, Xj and the following integral 



Uj c ^rU(X 1 ,...,X n )S(X 1 ,...,X n ) = (2.15) 



if the function S is symmetric for at least one pair of A-s. Therefore for an arbitrary function 
F n (Ai, . . . , A n ) one can transpose any pair of A-s taking into consideration appearance of additional 
sign because of the antisymmetry of U n (Xi, . . . , A n ). For example, if one transposes Xj with A& 
one gets 

F n {. . . , Xj, . . . , A fc , . . .) F n (. ...A,, X r .. .). (2.16) 

II. The reduction of the power of denominator for T n is based on two relations which can 
be checked directly 

111 



Xk — Xi — i Xj — Xi — i Xj — Xk — i 
1 1 + . 1 1 . 1 1 ^ 

Xj — Xi — i Xj — Xk — i Xk — Xi — i Xj — A; — i Xk — A/ — i Xj — Xk — i 
j'-i i j- 1 i j- 1 i 

^ \ j - \ k - i ^ i \.-\ k -i n Afc - Ai 

i jt k 
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In Appendices A and B we will show how the reduction can be performed for n = 3 and n — 4. 
Unfortunately, so far we have not succeeded in finding a result for general n. 

III. The ratio 

T n +i (Ai, . . . , A n+ i) _ II"=i(Aj + i) A" +1 ^ 
^n(Al, . . . , A n ) rij=i(An+i — Aj — i) 

is symmetric with respect to any permutation of Ai, . . . , A n . Therefore the relation (2.19) allows us 
to use the result T n also for derivation of T n+1 if this result was obtained by applying the relations 
(2.16-2.18) from I and II. 

IV. Proposition 1 

Let the function /(Ai, . . . , A n ) have only poles of the form 1/ (Aj — \ + id) with a an integer i.e. 
the product U n (Xi, . . . , X n )f(Xi, . . . , X n ) does not have poles of that kind. Then 

A}" /(-.., A,, . . .) ~ -(A, + i) m ./(.... A, +<,...) (2.20) 

where m is an integer and m>n. 

Proof 

Let us suppose that all variables Ai, . . . , A 3 -_i, Xj +1 , . . . , A n are fixed. Extracting from U n (Xi, . . . , A n ) 
the function which depends on Aj one gets 

Jc_ 1/2 2m sinh n 7rAj J J Jc_ m 2m sinh n 7rAj J J 

= -(/ "/ +/ j^ n^sinh^-A.) 

vc_3 /2 ic_ 1/2 Jc_ 1/2 2m smh 7rAj j j 

where w(Xj) = /(..., Aj, . . .). The first step here was to shift integration variable Aj — > Aj + i 
and to use the fact that sinh7r(x + i) = — sinli7rx. The two first integrals in the last expression 
are equal to a contour integral around the point Xj = — i in a complex plane of the variable Xj. 
Since, m > n the term (Aj + i) m which is in the numerator and corresponds to a zero of order m 
compensates the pole from the term sinh n 7rAj in the denominator. Therefore the contribution of 
those two integrals is zero and we immediately come to the statement (2.20). 

One can get two useful corollaries from the proposition 1. 
Corollary 1 

Xf g(X 1 , . . . , A,-, . . . , A n ) ~ i-i + ••••A, A„) (2.21) 

A k=0 

where the function <?(Ai, . . . , Xj, . . . , A n ) does not depend on Xj and as above it is implied that 
m > n. 

Proof 

The relation (2.21) is easy to derive using the relation (A™ - + (Xj + i) m )g(Xi, . . . , Xj, . . . , X n ) ~ 
or equivalently X 1 J l g(X 1 , . . . , Xj, . . . , X n ) ~ 1/2(A™ - (A, + i) m )g(X 1 , ....X, A„). 
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Corollary 2 

ym—l 

t ~ — r~ fi^Ai, • • • , Afc, . . • , Aj, . . . , A n ) ~ 

Ak — Aj 

( 171 \ jl\m-l 
m \ 7 J 1 A j m—1 

~ j— + E A i( A i + < ) m - 1 - , >( A i»---»^---»\v» A n) ( 2 - 22 ) 

where 

/ to \ to! 
W / Z! (to-/)! 

is binomial coefficient and the function g(\\, . . . , A*,, . . . , Aj, . . . , A n ) does not depend on \ k and Aj 
and to > n. 

Proof 

Using the proposition 1 we get 

A m ~ ~ (A ■ + i) m 
' <?(Ai, • • • , A*;, . . . , Aj, . . . , A n ) ~ — : <7(Ai, . . . , A fe , . . . , Aj, . . . , A n ) = 



A\ V *-L7 • • • 7 ' 'ft: 7 • • • 7 ' *J 7 • • • 7 ' *7t/ \ \ 

k ~~ A j A k ~ Aj ~ % 



= K, + )<?(Ai,---,A fc ,...,Aj,...,A n ) 



Afe — Aj — i Ak — Aj — 2 

j Xk + i) m | A^-(Aj + ^r- 
^ — Aj Afc — Aj 



(— : 7 1 — ~r r r — ) <?(Ai, . . . , Afc, . . . , Aj, . . . , A n ) 

V Ak, — Ao Ak — Aj — I ' 



or 



~ (-: : 1 : : — ) #(Ai, . . . , A fc , . . . , Aj, . . . , A n J ~ 

Ak — Aj Ak — Aj — I ' 

(A J+< r-Ar + Ai-(y <r )s(x x *„...ao~o. 



Xk — Aj Afc — Aj — i 

Then expanding both numerators according to the formulae 



in 

(A,+*r-A™ = £(7)^ 



=1 

m—l 



K ~ (A, + ^) m = (Xk - A, - E A ^ + *)"^ 

2=0 

we arrive at the formula (2.22). 

With the help of the corollaries 1 and 2 one can effectively reduce the power of the numerator 
in T n . 

V. For the calculation of integrals we need the following 
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Proposition 2 Let the integral fc_ 1/2 j^i el5X ^nh"wx ^ e convergent for any real number 
5 > and an integer N. Further, let the function /(A) be analytic in the whole complex plane A 
and satisfy the following two conditions 

(2.23) 
(2.24) 



lim \e- n * R f(iy - i/2 ± R)\ = 0, 

K— >oo 

Um , N f {x -i/2±iN) 

N^oo COsh 7TX 



where the first limit is uniform in y , when y G [0,7V]. The second limit is uniform in x for any 
real x. The value 8' is a fixed real positive number (5' > 0). Then 

dX /(A) 



,„ = lim S n) U) ^Ti-lf 1 e- 51 f(il + e) = 
2 2tu sinh n ttA 5^o+ °^ 

oo 

= - hm dW(e)_ 5:(-l) ,B e- a /(-« + e) 



z=i 



where a differential operator d( n )(e) looks as follows 



S n \e) 



n-l 



\ y( " 



n-l-i 



0<2fe<n 



(7TOJ 

(2fc + l)! 



In particular, for n = 2, 3, 4 



7i 2 oe 



Proof Let 



d< 3 >(e) 
F(X,N) 



1 ^ 

2 a.2i 



7T 2 

_L_ 9^ 

'3tt2 l 9e 4^2 9? 



2tt 
2 ,d 



)■ 



ciA 



,j<5A 



/(A) 



/Cjv 27ri sinh n 7rA 
where <5 > and Cat is a rectangular contour shown in Fig. 1. 



—Cn-i/2 , 



(2.25) 
(2.26) 

(2.27) 

(2.28) 
(2.29) 
(2.30) 

(2.31) 



, (JV-l)i 



C: 



(2 



C-l/2 



Fig. 1 The contour Cat 
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The contours C_i/2 and — Cat-i/2 correspond to lower and upper horizontal parts of the contour 
Cn respectively (sign — is because the contour Cjv-1/2 should be taken in the opposite direction). 
The contours C+ and correspond to the right and left vertical parts of the contour and 
have real parts +00 and —00 respectively. Due to the Cauchy theorem one has 



JV-l 



F(S,N) = F_ 1/2 - F N _ 1/2 + F + + F_ = rf (n) (e)_ ^ (-l) ln e~ 5l+i5e f(il + e) (2.32) 



1=0 

where 



Fn-1/2 



= f dA cl5x J (A) 

Jc_ 1/2 2iri sinh n 7rA 
= f d\ cl5x /(A) 
Jc N _ 1/2 27ri sinh™ 7rA 

^ = 4^& (2 - 35 » 

The r.h.s. of the formula (2.32) is a result of the calculation of residues corresponding to zeros of 
the denominator sinh n 7rA which are placed inside the contour Cn- 
The first step is to prove that the integrals over the contours 

F± = 0. (2.36) 

Actually one has 

d\ iSX /(A) 



F ± = lim / — e iS 



sinh n 7rA 



- li£L±C HF^W^ e " ±R+ '"' /2,/(±fl + iy - i/2) (2 - 37) 

where vertical contours C\ are defined as follows: {A = ±R + iy — i/2; y G [0, A 7 "]}. Then 

|F ± | < — lim ( N — $ e^l/fitf + - i/2)| = 

1 ±! - 2tiR-ooJ sinh^Tri? V y ' n 

= B lim / n ^n ^i/^ + *f - */2)| (2.38) 

where we have used a simple fact that for R > 

I cosh7r(±i? + > sinh7nR. 

The uniform character of the limit (2.23) allows to interchange the order of the integration over 
y and the limiting procedure R — > 00. Indeed, the condition (2.23) means that for any small real 
number e there exists a real R e which is independent of y such that for any R> R e 

e- n7rR \f(±R + iy - i/2)\ < e. 

Therefore for R > R e 

r N dnp~ s y p r N 

/ , o ^ e- n7rR \f(±R + iy - 2/2)1 < — — / dye~ 5y = 

Jo (1 - e- 2 * R ) n UK y ' n (1 - e~^ R ) n Jo y 
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e -N8 e l _ e -N8 

< 



Hence we have got that 

r N dye~ 5y 



lim / . y€ - PN e-™ R \f(±R + iy - i/2)\ = 



and we come to the statement (2.36). 

The next our step is to prove that for a fixed real 5 > 



lim F N _ 1/2 = (2.39) 



Indeed, 

iS(\ + iN) 

f(\ + iN). 

N ~ 1 I 2 Jc N _ 1/2 2ni sinh n 7rA Jc_ 1/2 2m sinh n 7rA 



x j 2 Olllll /\ J ^-l/2 

e- SN , r dX e- 5iV /■«> f/.c 



2tt 



/ -T^T /( A + ^ V- / — S— - </2 + <JV)| (2-40) 

Jc_ 1/2 smh" ttA v 71 2tt J-oo cosh" ttx 1 v ' 71 v y 

As above due to the uniform character of the limit (2.24) one can interchange the integration over 
x and the limit iV — > oo. Therefore the last expression in (2.40) tends to zero when N — > oo and 
we come to the statement (2.39). 

So for a fixed 5 > we conclude from (2.32) that 



lim F(6,N) = [ ^e iSX - = S n \e)^ V(-l) fa e~ sl+iSe f(il + e). 

iv^oo v ' 1 Jc_ ll2 2m smh n 7rA 



-1/2 



The convergence of the sum in r.h.s. of the last expression is guaranteed by the convergence of the 
integral in the l.h.s. Finally, after taking the limit lim,5^ + we come to the formula (2.25). Let us 
note that generally speaking we can not interchange the order of the limits — > oo and 5 — > + . 

Another form (2.26) can be proved in a similar way if considering a contour which is analogous 
to Cn but placed in a lower half-plane and a real number 5 should be taken negative. It completes 
the proof of the proposition 2. 



Using the proposition 2 we can get the following 
Proposition 3 

II / irrUn(Xi, . . . , A„) F„(Ai, . . . , A n ) = F n (ei, . . . , e n ) (2.41) 

j = l ^1 

where the multiple integral is convergent and the product U n (Xi, . . . , X n ) F„(Ai, . . . , A„) does 
not have any other poles besides the poles of the denominator lTj=i sinh n nXj of the function 
U n (Xi, . . . , A n ). The function 

G^(A 1? ...,A„) - n^.^n^ F n (Ai, . . . , A n ) (2.42) 



11 



should satisfy conditions which generalize (2.23) and (2.24) 

lim \e- n * R G®(\ 1 ,...,\ j - 1 ,ix j -i/2±R,\ j+1 ,...,\ n )\ = 0, j = l,...,n (2.43) 

lim \e~ SN G^\\i, . . . , Xj — i/2 ± iN, Xj+i, ■ ■ ■ , A n )|/cosh n irxj = 0, j = l,...,n (2.44) 

N — >oo 

where for each j both limits are uniform on real numbers xi, . . . , x n with A& = Xk + im k /2, k ^ j 
and some integers m&. For the first limit Xj G [0,N] while for the second one Xj is any real number. 

D (n) 

is a differential operator 

^")=7r^n^ (e .)^ o JJ sinh^-e,) (2.45) 

j=l l<k<j<n 

and 

oo oo 

F n ( ei ,...,e n ) = lim V)(_i)'i e -*i'i.. hm ^(-i^e-^F^ih + e u . . . ,il n + e n ) (2.46) 

5l ^ 0+ l 1= 5n ^ l n =0 

where each sum X^°° =0 (— l^e^^F^Ai, . . . , Xj-\,i lj + €j, Aj+i, . . . , A n ) should be convergent uni- 
formly in other arguments Ai, . . . , Xj-i, Xj + \, . . . ,X n (Xk = x k + im.k/2, k ^ j). 

Proof The formula (2.41) can be got by the recursive application of the formula (2.25) to 
each integral in the l.h.s. of (2.41) and by taking into account the manifest form (2.3) of the 
function U n (Xi, . . . , X n ). After j — 1-th application one can interchange the integration over Xj 
with the sums over h, . . . , lj-i due to the uniform convergence of these sums. Due to the conditions 
(2.43) and (2.44) one can apply the proposition 2 when integrating over the variable Xj. 

Let us note that each sum X^ =0 ( — e~ s i lj F n (. . . ,iXj + €j, . . .) in the formula (2.46) can be 
substituted by — J2i°=i(~ e~ Sjlj F n (. . . , — iXj + €j, . . .) corresponding to the choice of the contour 
in the lower half-plane. In what follows we will use this fact depending on a convenience. 

Let us consider a special class of functions F n (ei, . . . , e n ) that does not have a singularity when 
6j — > for j = 1, ... ,n. In this case one can expand F n (ei, . . . , e n ) into the infinite series on 
powers of e-s. We have checked that for n < 4 the differential operator given by (2.45) 
when acting on some monomial e^ 1 . . . e^ 1 makes non-zero contribution only for monomial of the 
form e° (1)6^(2) • • • € ™7n) where a is some element of the permutation group S n of n elements. More 
precisely, we can write 

D *" = Trh n (s-s-) (2 - 47) 

1 ij=l J ! l<k<j<n V ° e k 6e 3 ' 6-0 

We proved it for n < 4 and we use this formula only in this case. We believe that this relation is 
valid for any n but this fact is still to be proven. 

VI. Now let us discuss the integrals of a special form, namely, when the function F n (Ai, . . . , A n ) 
from the proposition 3 is rational on it's arguments and the function 

Y[ sinh7r(A i - A fc )F n (Ai, . . . , A n ) 

l<k<j<n 

is analytic. It means that the function F n (Xi, . . . , A n ) can have only simple poles when Xj — > 
Xk + mi with an integer m. As it is seen from the definitions (2.4), (2.7) both the function 
T n (Ai, . . . , A n ) and the function T^(Ai, . . . , A n ) are of that form. 
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Let us show that the conditions of applicability of the propositions 2 and 3 are fulfilled for such 
a function F n (Xi, . . . , A n ). To show that the first condition (2.43) is satisfied let us fix without loss 
of generality j '• = 1 and take into account that < x\ < N. First let us consider the case when 
the rational function F n does not have poles at all. Then if m 2 , . . . , m n have the same parity 

e~ nwR \G<P(ix! - i/2 ± R,x 2 + im 2 /2,...,x n + im n /2)\ = 



n 



= e 



—nwR 



j [ | sinh 7r(x k =F R — ixi 



| IIi<fc<Kn I sinh7r(rg; - x k ) 

rifc=2 COS 11 " ™k 



k=2 llfc=2 

\F n (ixi - i/2 ± R, x 2 + im 2 /2, . . . ,x n + im n /2)\ < 
< e -n,n 2 i^K^ ft I sinhvr^J R - ^x l )\ _ ±R ^ + ^ + < 

k=2 COSh 7lX k 

< e _ BTfl 2 (n^H ft cosh7r(x fc =F g) _ . /2 ± fl>a . a + , m2/2; . . . >a . w + imB / 2 )| < 



fc=2 



cosh 7ra; fc 



D (n-2)(n-l) , JL 1 

< e - n7ri? 2 2 \ cos\mR + sinhTri?)"" 1 J] 



^ COSh 7TX k 

F n (ixi - V 2 ± R,x 2 + im 2 /2, ...,x n + im„/2)| = 



(n-2)(n-l) „ 

2 5 e" wfl |F n (ia:i - V 2 ± R,x 2 + im 2 /2,...,x n + im n /2)\ (2.48) 



where we have used an inequality 



Y[ | sinh7r(x/ - x k ) 

ITl-^MTl-ll -a r j 

< 2 — ~ II Z3^-T (2-49) 



72 

Kfe<Kn ^ (n-2)(n-l) 



n C Osh ra 7TX fe fc=2 COSh TTX* 



which can be checked directly. For a set of arbitrary integers m k , k — 2, . . . , n one can repeat the 
derivation above also. 

Since < X\ < N and F n is rational without poles one can find real numbers R* and M(N) > 
which is independent of R such that for R > R* 

n y 

n — : |F n (ia;i - i/2 ± R,x 2 + im 2 /2,...,x n + im n /2)\ < R S M(N) 

^ = A 2 cosh -nx k 

with some power s. Hence, we get that the expression (2.48) can not exceed 

(n-2)(n-l) D 

2- — 2 — 1 e ~ nR R S M(N) 

which tends to zero when R — > oo independently of the variables xi, . . . , x n and we get the uniform 
character of the limit (2.43) on these variables. 

Let us prove the validity of the second limit (2.44). Indeed, for a given integer j — 1, . . . n and 
X k = x k + im k /2, k 7^ j with integers m k of the same parity the function 

\e~ SN GV>{\ u ...,\ j - 1 ,x j - t/2 ± iN,X j+1 ,...,X n )\/cosh n 7rx j 
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is bounded for any real numbers x k € (—00, 00) and a positive integer N. Therefore one can use 
again the inequality (2.49) in order to get 

tt 1 • v. / .. n |ni<ifc<KJ sinh7r (^ Ili<jfc<Kn I sinh 7r(rr t - x k )\ 

nismh^-^-^i — ^_ . rot — = - is^ih^ — £ 



fc=2 



n(n-l) I 
< 1—^- 



Ilfc=l cos h TT^Jfc 

Again one can repeat this for arbitrary integers m 2 , . . . , m n . 

Since F n is rational the maximum over the real variables xi, . . . ,x n of the function 

\F n (xi - i/2 + iN,x 2 + im 2 /2, . . . ,x n + im n /2)\ 



Ilfc=l cosn TT^fc 

can be some power of N, say, N s ' multiplied with some constant which is independent of iV when 
N > N* with N* is a big enough integer. Therefore we get the limit 

lim e- 5N N s ' = 

N— >oo 

for any real 5 > independently of variables x±, . . . ,x n and we come to the uniform limit (2.44). 

Suppose the function F n (Xi, . . . , A n ) has a simple pole of a type 1/ (A& — A; + iafcj) with an 
integer clm- Let us restrict ourself only with the case a k i = because only such poles can appear 
in the expression for a canonical form (2.7). Then we can write 



j Sinh7r(A fc - A;) < 

Afc — A; 



it sinh 1 if I Afe — A; I < — 

T 

7r| sinli7r(A fc - A*)| if |A fe - X t \ > - 

71 



and use again the technique described above. 

Let us comment on a question of uniform convergence of a sum 

00 

X] (-l) lie ~ Sjlj Fn(^i, Aj-i, i lj + e j ,X j+1 ,...,X n ) 

as well as how to proceed further by considering two typical examples. Let us take for simplicity 
n — 2. The integral we need looks as follows 

W £?/ ^ S ^t 2 ^\ (2-50) 
Jc_ 1/2 Atii Jc_ 1/2 2,711 sinh 7rAi smh 7rA 2 

i As a first example let us consider again the case when F 2 does not have poles at all i.e. 
F 2 (Ai, A 2 ) is some polynomial on Ai and A 2 . First let us integrate over Ai using the formula (2.25) 

J 2 = / , , 2 lim ^ + ^(6 1 ) £l ^ smh7r(A 2 -6 1 )E(- 1 )' 1 ^ 1 ^2(^ + 6 1 ,A 2 ) (2.51) 

JC_ 1/2 2m smh 7rA 2 5i->o+ 1 fi^ Q 

Since F 2 (Ai, A 2 ) is a polynomial then F 2 (i/i + ei, A 2 ) is a polynomial on Zi as well. Let us pick out 
some monomial on l x from it, say, 

Z?Fx(A 2 ) (2.52) 
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where a is a non-negative integer and -F\(A 2 ) is a polynomial on A2. Actually it has a factorized 
form. Therefore in this case we do not have any problems with the uniform convergence and the 
corresponding contribution into J 2 is as follows 

/ 9 , , , 2 2 Fi(A 2 ) lim S 2 \e l ) e ^ smhn(\ 2 -e l )Y,(-l) h e- Sll H a 1 = 
Jc_ 1/2 2m smh 7rA 2 <5i->o+ 1 £r Q 

d\ 2 



= L„ 5— ^^(^)^( ei )«_„Bi n l-(A 2 - « lW .) (2.53) 



where 



00 



p(a) = lim VK^IJ (2.54) 



Let us adduce a number of first values of p(a) 

P(0) = ^ AK1) 1, P(2) = 0, aK3) = |, P(4) = 0, p(5) = ~ (2.55) 

Since Fi(A 2 ) is a polynomial we can treat the integral in (2.53) in a similar way as the integral 
over Ai. In the very end we should calculate the limits e\ — > and e 2 — > and get the final answer 
for J 2 . 

m The second example corresponds to an existing of a simple pole, namely, when 

F 2 (Ai,A 2 ) = - g( ^ Aa) (2.56) 
A 2 — \\ — ia\2 

where Q(\i, A 2 ) is a polynomial on Ai, A 2 . As above we shall consider only the case ai 2 = 0. So 
in this case doing the first integration one gets 

h = I 9 % X ,lim dW( f0a ^ a ^ T(A ,- fl )f : (-iy. e -».i. q("- + t -^) (2.57) 
ic_ 1/2 27n smh 7rA 2 <5i— >o+ 1 A 2 — ill — e 1 

Since Q{il\ + ei, A 2 ) is a polynomial on Zi also again let us pick out some monomial from it 

with an integer a' > and a polynomial Q'(A 2 ). Then the corresponding contribution into the 
expression (2.57) for J 2 looks as follows 



r rl A 00 l a ' 

/ -— g'(A 2 ) lim ^(ex^sinhTrCAa-ex) £ (-l)^ 1 — -i (2.58) 

ic_ 1/2 2m smh 7rA 2 ii^o+ 1 ^ A 2 — z/i — ei 

Since A 2 = x 2 — i/2 with a real number x 2 G (—00, 00) the denominator 

1 i 



Xi—ih — e-i h + 1/2 + ix 2 — ie\ 



and Re(li + 1/2 + ix 2 — iei) > 1/2 because l\ > 0. Therefore one can use an evident integral 
representation 

l -= i l *V (2.59) 
a Jo s 
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which is valid if Re(a) > 0. Hence the sum in the expression (2.58) is 

"i ds 



' ' fl W Q 

i y;(-i) ,i e-* i,i zf / ^£^+i/2+« 2 - iei (2 60) 



Zi=0 



Since for Si > the sum 

where JV' is a positive integer converges uniformly in s G [0, 1] when N' — > oo then one can 
interchange the sum over Zi and the integration over s. The result for (2.60) is as follows 

r 1 r) , 1 /■! a ' +1 1 

i / ™ 8 »*-*i(°f 1 = i f ^-fci ^ C(a',fc)- ( 2 - 61 ) 

7o s <9<V 1 + e~ Sl s Jo s j£l (! + e- 5l s) fc 

where C(a', fc) are some rational coefficients. So the contribution of /c-th term into the expression 
(2.58) is 



'c_ 1/2 2m sinh 2 7rA 2 
It is not very difficult to check that the function 

Q / (A 2 )sinh7r(A 2 -e 1 )s lA2 
satisfies the conditions (2.23) and (2.24) of the proposition 2 with n = 2. Moreover the integral 



/ ^e^ Sinh7r( ^ 2 ~ ei) g-(A 2 )^ 
Jc_ 1/2 2m sinh 7rA 2 



1/2 

converges uniformly in s £ [0,1] . Since it, actually, does not depend on Si also and S 2 \e\) 
given by (2.28) is the second derivative on e x one can calculate the integral over A 2 first using the 
formulae (2.25-2.26) of the proposition 2 and then calculate the integral over s and take the limits 
on ei and Si. 

Therefore we have shown that the formula (2.41) is correct for n — 2. Below we will use this 
formula for the concrete calculation of P{2). 

A generalization of our discussion to the arbitrary n is straightforward and we will not do it 
here. 

The most efficient way of taking the integrals is as follows. Fist we apply the formula (2.41) 
using either (2.25) or (2.26) in such a way that the denominators like 

becomes 



A 2 - Ai h + l 2 + i(e 2 - €i) 

with li > and l 2 > 1. In this case we shall not face a singularities like l/(e 2 — ei) and the whole 
expression will be analytic on ei,...,e n . Hence, we can use the differential operator (2.47). Then 
after using the formula like (2.59) one can get rid of all such denominators and expand over e-s 
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until the order which still makes a non-zero contribution after applying the differential operator 
(2.47). It can be applied after taking all summations in (2.46). As the last step one should take 
integrals on auxiliary variables like s appeared above. 

Below we shall also take all 5i,...,5 n to be zero at once implying the limiting procedure 
lim s ._> + described above. 

Now let us illustrate how the whole procedure works for a simple case P{2) 

P(2) = tt 3 [ p:l p: Aa) (2 . 63 ) 

Jc_i/2 27ri jc*_ 1/2 2m sinh 7rAi sinh 11X2 

In this case it is very simple to perform the first step, namely, to get the "canonical" form (2.7) 
described in the beginning of the Section because we do not need to reduce a power of denominator 
in this case. Indeed, 

r ! (A 1 ,A 2 ) = ^i±^ = A 1 + ! + i^±^~A 1 --ii- (2.64) 

A2 — Ai — l A2 — Ai — l A2 — \\ 

where we have used the property I and the formula (2.20) from the proposition 1 of the item IV. 
Then using the formula (2.22) of the corollary 2 for m = 2, 3 one gets 

+ X\ +X 2 (X 1 + i) + (Ai + z) 



A2 — Ai 3 V A2 — Ai 

iAi 1 1 i,., ,11 11 1,1,11 
- 77(^1) ~ ; 77^ r- + 77 Ai = -Ai + 



A2 — Ai 3 A2 — Ai 3 2 A2 — Ai 3 A2 — Ai 3 3 6 A2 — Ai 

Substituting it to the formula (2.64) we get 

T 2 (A!,A 2 ) ~ T 2 C (A 1? A 2 ) (2.65) 

where 

T 2 c (A l5 A2) = Ui + Iy^V ( 2 - 66 ) 
3 6 A2 — Ai 

and this is the "canonical" form for T 2 i.e. the polynomials P and P± from (2.7-2.8) are equal 
to 4/3Ai and 1/6 respectively. 

Let us take the integral from the first term using the formula (2.46), our comments about the 
limiting procedure like in the formulae (2.25-2.26) of the proposition 2 

4 2) = - 3 I d -±[ ^ sinh^- AQ 4 = (2) 4 g , g h = 

Jc_ 1/2 2m Jc_ 1/2 2m sinh 2 vrAi sinh 2 ttA 2 3 3 ^ ' v ; ^ ' 

= B «^ M1)+flp( o )W o, = o^-L + ^ = ( ^-^J ( -i + |) = i (, 67) 

where we have used the formulae (2.55) for p(b) given by (2.54) implying the limiting procedure 
as it was explained above. 
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The second integral is treated as it was described in item VI with the help of the integral 
representation (2.59) 



(2) 3 f dX 1 f d\ 2 sinh7r(A 2 - Ai) 1 1 



c_ 1/2 27ri Jc_ 1/2 27ri sinh irXi sinh 7rA 2 6 A 2 — Ai 

(_-t\ OO OO / -1^1+^2 • OO OO / l\il+I 2 

D (2) (_iZ V V ' } = D {2) - V V 



6 z 1= i ; 2=0 ^2 + iZi + e 2 - ei 6 ^ h=Q h + l 2 + i(e t - e 2 ) 

= D (2) % _ f^h+hsh+h+i^-e,) = D (2) L^l ds * = 

6^0 S ^ 6 7o (1 + S ) 2 

. 9 9 , H) , s^i"^) 1 /-i Ins In 2 

= fe-^^o^X ds 0TW = 3io rfs (TT^ = — (2 ' 68) 

Summing up two answers (2.67) and (2.68) we get the result 

P(2) = Jf + Jf = 5 - (2.69) 

which coincides with the formula (1.7). 

In the Appendices A and B we shall derive the formulae (1.8) and (1.9) for P(3) and P(4) 
respectively. In the end of this Section let us note that both results (1.8) and (1.9) are expressed 
in terms of the logarithmic function and the Riemann zeta function of odd arguments and do not 
depend on polylogarithms like, for example, Li 4 (l/2). All coefficients before those functions in 
(1.6-1.9) are rational. Also they do not contain any powers of n which could be considered as 
Riemann zeta functions of even arguments, see the formula (1.4) from the Introduction. 

Our conjecture is that the final answer for any P(n) will also be expressed in 
terms of logarithm In 2 and Riemann zeta functions ((k) with odd integers k and with 
rational coefficients. 

In fact, this conjecture is intimately connected with our hypothesis from the beginning of 
this Section that the function T n (2.4) can be reduced to the "canonical" form. Looking at the 
"canonical" form (2.7) one can conclude that only Riemann zeta functions and their products can 
enter into the final answer because all the denominators in the r.h.s. of (2.7) are split out. It 
means that after applying the formula (2.46) the multiple summation can be performed by pairs, 
say, Z)j 2fc _ 1 and J2i 2k - Each pair of these summations results in some combination of zeta functions. 



3. Conclusion 

We want to emphasize an interesting connection between integrable and disordered models. In 
order to describe correlations in integrable models one can use integrable integral operators [23]. 
On the other hand Tracy and Widom showed that in matrix models the distribution of eigenvalues 
and level spacing can be described by the integral operators, belonging to the same integrable class 
[22]. 

Our current work supports this link between integrable models and chaotic models. Riemann 
zeta function appears in the description of both kind of models. 
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Let us repeat that the main result of this paper is the calculation of P(3) and P(4) (1.8- 
1.9) by means of the multi-integral representation (2.1). The fact that only the logarithm In 2 
and Riemann zeta function with odd arguments participate in the answers for P(l), . . . , P(4) 
and with rational coefficients before these functions allows us to suppose that this is the general 
property of P(n). One could compare the calculation of P(n) with the many-loop calculation of 
the self-energy diagrams in the renormalizable quantum field theory which can also be expressed 
in terms of ( functions of odd arguments [1] . 

Unfortunately, so far we have not got even a conjecture for P(n) but we believe that it is not 
an unsolvable problem. May be already after calculation of P(5) one could guess the right formula 
for a generic case P{n). It would give an answer to the question discussed in the previous section, 
namely, the question about the law of decay of P(n) when n tends to infinity. 

Also it would be interesting to generalize above results to the XXZ spin chain. Some interesting 
conjectures were recently invented by Razumov and Stroganov [13] for the special case of the XXZ 
model with A = —1/2. These conjectures would be supported if it were possible to get P(n) from 
the general integral representation obtained by the RIMS group [7]. 
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5. Appendix A 

Here we discuss in detail the calculation of P(3) performed by means of the general procedure 
described above. 

As was pointed out in the beginning of Section 2 the first step should be a reduction of the 
function T^(\i, A 2 , A3) to the form (2.7) which we have called the "canonical" form. In comparison 
with the case P(2) here we should reduce the power of the denominator in T3(Ai, A2, A3). To do 
this we will use the formula (2.17) from the item II. Namely, 

T (\ \ \\ _ (Ai + i) 2 (A 2 + i)A 2 A^ (3) (3) (3) 

(A 2 - Ai - i)(\ 3 - Ai - i)(\ 3 - A 2 - i) 



where 



,(3) _ (A 1 +i) 2 (A 2 + i)A 2 A| 

1 " ? (A 3 -A 1 -0(A3-A 2 -0' [ } 



7 (3) . . . (Ai + i) 2 (A 2 + i)A 2 A 2 

(A 2 - Ai - i)(A 3 - Ai - i)' 
7 ( 3 ) . : (Ai + 2) 2 (A 2 + QA 2 A 2 

(A 2 - Ai - i)(\ 3 - A 2 - i)' 



(5.3) 
(5.4) 



(3) 

Due to the 1 <-> 2 symmetry of the denominator the first term I\ can be simplified as follows 

(3) (Ai + i)(A 2 + i)A 2 A 2 (Ai + i)(A 2 + i)A 2 

/ — ■ 



(A 3 - Ai - i)(X 3 - A 2 - i) A 3 - Ai 
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= (Ai + 2)(A 2 + 2)(Ai + A 3 + «) + ^ — J -f — ~ A 2 A 2 - ^- — \ K . (5.5) 

A 3 — Ai — i A 2 — Ai — % 

The denominator of the second term (5.3) has the symmetry under the transposition 2 <-> 3. 
Therefore it can also be simplified 

r(3 ) (Ax + i) 2 A 2 A2 (A 1+ i) 2 A 2 A 3 

(A 2 - Ai - i)(\ 3 - Ai - i) A 2 -Ai-i 

<\ (Ai + Q 3 A 3 X 2 X (Ai + ^) 3 A 3 

= -(Ai+i) A 3 - : ~ A : A 2 - : (5.6) 

A 2 — Ax — l A 2 — A\ — % 

The third term (5.4) is treated as follows 

if = -i(X l + t) 2 (X 2 + t)(X 3 + X 2 + t) - 
^.(A 1+ *) 2 (A 2 + *) 3 . (A 1+ *) 3 A 2 . (A 1 + *) 3 A 3 



A 3 - A 2 - i X 2 — Xi—i (A 2 - Ai - i)(A 3 - A 2 - i) 

\2\ (A 1+ Q 3 (A 3 + .) 2 (A 1+ Q 3 A§ (Ai+*) 3 Aj 
~ " 1 2 X 2 -X 1 -i + ~ t (A 2 -A 1 -*)(A 3 -A 2 -») (5 ' 7) 

Now adding up all the three terms together we get 

t<\ \ \\ \*\ (Ai+^) 3 Aj x2 , (Ai + 3 A 3 

J 3(^1, A 2 , A 3 J ~ A X A 2 — — 7T : rr ~ — A 2 A 3 



(A 2 - Ai - i)(A 3 -X 2 -i) (A 2 - X 1 - i)(X 3 - A 2 - i) ' 

(5.8) 

Let us note that up to this moment we have used only the symmetry property (2.16) from the 
item I, the formula (2.17) from the item II and a simple algebra. 

Now we would like to use the formula (2.20) of the proposition 1 for m = 3 and again apply 
the transposition formula (2.16) 



^~3(Ai,A 2 ,A 3 ) ~ — A 2 A 3 



■> . a 3 (A ;! • iy 



(A 2 _A 1 )(A 3 -A 2 ) 

, , . (A x +i) 3 Al 2 .A?A^ + 3iA 2 A|-3AiA|-iA^ 

" t (A 2 -A 1 )(A 3 -A 2 ) " + * (A 2 -A 1 )(A 3 -A 2 ) 



2 3A 2 A 2 + 3iAiA 3 (A 3 + A 2 ) — A 2 — A 3 A 2 — A 2 , 
— a 2 a 3 — 



A 2 — Ai 

2 3A 2 A 2 + MAiA§ + 3*A 2 A 3 - A 2 - A3A1 - X\ 
~ -A 2 A 3 (5.9) 

A 2 — Ai 

Now we can reduce the power of Ai in the numerator of the second term (5.9) by applying the 
formula (2.22) of the corollary 2 from the item IV. Doing this we finally get the "canonical form" 
(2.7-2.8) of T 3 

p(3) 

T 3 (A 1? A 2 ,A 3 ) ~ T^Ax^As) = P (3) + ( 5 - 10 ) 

\ 2 — \\ 
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where the polynomials Pq 3 ^ and Pi 3 "* are as follows 

P W = -2\ 2 \ 2 3 , p{ 3) = I - iAx - i\ 3 - 2AxA 3 (5.11) 

Let us note that if we express variables Xj through the real variables Xj via Xj = Xj — i/2 in order 
to get the polynomials Pf^ ( see the formula (2.11)) we get especially simple formulae, namely, 



So, the function 



P {3) = -2x 2 xl p] 3) = ~- ^x 3 . (5-12) 

6(3) 

T^ Xl ,x 2 ,x 3 ) = P (3) + (5.13) 

x 2 - Xi 

is odd i.e. 

f 3 c (-x 1 ,-x 2 ,-x 3 ) = -f^(x 1 ,x 2 ,x 3 ) 

as it should be according to the formula (2.13) from the beginning of Section 2. 
Now we are ready to calculate the integral in order to get the result for P(3) 

3 „ j\ 

p (3) = II / T 4tf3(Ai,A 2 ,As)r 3 (Ai,A 2 ,A3) = 4" ] + J? ] (5-14) 

j=l "^C-l/2 

where 

4" = Uj C m £ *<A, A, A 3 ) if >, 7f» = n/ c i/2 i WA, A, A 3 ) ^ (5,5) 
Using the formulae (5.11), (2.46), (2.47), (2.55) we get 

= f[ / ^^(Ai,A 2 ,As)(-2)A 2 Al = 



_Y * C— 1/2 27TZ 

oo oo oo 

h=o 2 2 =0 « 3 =o 



oo oo oo 1 H f) 1 

= ^E(-i) h E(-i)' 2 E(-i)' 3 (-2)^ = I n (^-^) (-2)^ = ^ 

i 1= Z 2 =0 Z 3 =0 ° l<fc<j<3 atfe Ut 3 e->0 ^ 

(5.16) 

To calculate the second term jj 3 - 1 we should also use the integral representation (2.59) 



J=i J C-i/2 2vn A 2 - Ai 



L) (3)g ( 1)h g ( ~ »H*i + O - <(«3 + cs) - 2(-tZ 1 + Cl )(iZ S + C3) 



, , ^ , . (j il 2 + ih + e 2 -e 1 
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= ipWf^j If f]( l)faf)( 1) ^ ~ !l + k ~ lei ~ ^ ~ 2(Zi + ie 1 )(Z 3 -ie 3 ) 



h=i i 2 =o i 3 =o 



h + h + - £2) 
= <JW g (-!)■■ £ (-!)■■ -H 7^' +' lfc » 7 ^ = 

ii=l «2=0 <1 +<2 +H e l - £ 2j 

7o s 2 1= i/ 2 =o 12 



,(l + s) 2V 12 0/ ° y (l + s) 

(9 9 \ r 1 ds { 1 „ „ o 9x / .x In 3 s 



^3 



2 /0* a \ f 1 ds ( 1 o 2 m s , , N 

o ("» "a ) / 7-— n 2 (""" (- 3e i e 2 + 3eie 2 ) (-2) -— + 1 Ins 

2 l<fc<i<3 V ' 9 ^ d€ 3 J e-^0 J ° (1 + S) 2 V12 6 

= r'ds-^l^ _ 1 f'dsJ^I^ = _i n2 + ^(3) (5.17) 
7o l + s 2 12 Jo l + s 2 8 SV ; v ; 



1^ 

/o ~"(l + s) 2 12 Jo ""(l + s) 
Summing up Jo 3 " 1 and j[ 3 ^ we get the final answer (1.8) 



P(3) = J< 3 > + J< s) = i - In 2 + |C(3). (5.18) 

6. Appendix B 

Here we will calculate P(4). This case is much more complicated than the previous ones. But 
using our general technique we will try to simplify our discussion as much as possible. 

As above the first step of our scheme is to get the " canonical" form for the function 
T 4 (Ai, A2, A3, A4). 

But before start doing this let us list some useful formulae which can be derived from the 
proposition 1 and corollary 2 in the case n = 4 



A 3 -fa 2 + Ax + \ 1 



+ 0^2 + 7-A2 (6.1) 



A2 — Ax A2 — Ax 2 2 

Af — A? — iAx + 77; l,,o 2, . . 

— ^ + -A 3 + -AxA 2 + -A, (6.2) 



A2 — Ax A2 — Ax 5 5 5 

% .. -T - T2 , K x3 4. x3 , 13 x2 , 4. A I 



+ -A1A5 i\l + — A 2 + -i\ 2 (6.3) 



A2 — Ax A2 — Ax 3 3 6 3 3 

^i . -|A 2 - I^Ax - I 



A2 — Ax A2 — Ax 



+ 



+ ^A 3 - ^AxA 3 + ^A 3 + ^AxA 2 - ^A 2 + ^A 2 (6.4) 

(Ai + *) 3 „ M - 2Ax - |i + 1 A 2 + (6 5) 

A2 — Ax A2 — Ax 2 2 
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. - A? .-"- +A + + i Al A> + 2^ - h 2 (6.6) 

A2 — Ai A2 — Ai 



A 2 - Ai A, -A, +3^+3^ + «^--8^-3^-3 < 6 - 7 > 

Here we have omitted an arbitrary function g>(A 3 , A 4 ) that can be multiplied on the r.h.s. and 
l.h.s. of (6.1-6.7) as in the formula (2.22). 

Now we can start our derivation. Fortunately, due to the formula (2.19) of the observation III 
which in this case looks as follows 

T 4 (A 1 ,A 2 ,A3,A 4 ) _ (A! + i)(A 2 + i)(A 3 + i)A| 

T 3 (A 1; A 2 ,A 3 ) " (A4-A 1 -i)(A 4 -A 2 -i)(A4-A 3 -i) l ' J 

and has a symmetry under any permutation of the variables Ai, A 2 , A 3 we can use the result (5.8) 
that was obtained by means of the symmetry and simple algebra. So, after the application of the 
formulae (5.8) and (6.8) we can start with the following expression 

Tn x x xA (Ai + i)(A 2 + i)A 2 (A 3 + i)A§At 

J4(Ai, A 2 , A 3 , A 4 J ~- 



(A 4 - Ai - i)(X 4 - A 2 - i)(X 4 - A 3 - i) 



_ • (Ai+Q 4 (A 2 + t)(A 3 + i)A|Ag 

l (A 2 -Ai-i)(A 3 -A 2 -i)(A4-Ai-i)(A4-A 2 -i)(A4-A 3 -i) 1 ' ' 

The denominator of the first term in the r.h.s. of (6.9) is symmetric under permutation of Ai, A 2 , A 3 . 
Therefore it can be simplified as follows 

(Ai + i)(A 2 + i)A 2 (A 3 + i)A|Al (Ai+i)(A 2 + i)A 2 (A 3 + i)(A 3 + A 4 + i)A| 



(A 4 — Ai — «)(A 4 — A 2 — i)(\4 — A 3 - i) (A 4 — Ai — i)(A 4 - A 2 - i) 

(Ax + i)(A 2 + i)(A 3 + i)(\ 3 + A 4 + (Ai + i)(A 2 + i)(A 3 + i)\ 3 \\ 



A 4 — Ai — i A 4 — Ai — i 

(\ j_i\<\ \3 (A 2 + 2) (A 3 + i)A 3 A 4 2 3 A 4 (A 3 + i)(A 4 + i)A 4 

= (A 2 + ^J(A 3 + ?)A 3 A 4 ; ~ A 2 A 3 A 4 — — ~ 

A 4 — Ai — 1 A 2 — Ai + % 

2 3 Af(A 3 Al + ^-A 4 ) (A 1 +Q 4 (A 3 Al + al-A 4 ) 

~ A 2 A 3 A 4 - — — ~ A2A 3 A 4 H ~ 

A2 — Ai + 1 A 2 — Ai 

8, , 2X 3 (Aj + iAi - ^) (A 3 A1 + iAl - A 4 ) 

~ -A 2 A 3 A 4 - ^ _ ^ (b.lU) 

The latter formula was obtained with the help of the formula (6.6). In fact, it is nothing but the 
"canonical" form (2.7) for the first term in (6.9). 

Now let us treat the second term. Using the formulae (2.17) and (2.18) we can write it down 
as follows 



—1 



(Ai + 4 (A 2 + 0(A 3 + 2)A|Ai 



(A 2 - Ai - i)(X 3 - A 2 - i)(A 4 - Ai - i)(A 4 - A 2 - i)(A 4 -X 3 -i) 
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11 11 



2 (A 3 - A 2 - i)(A 4 - Ai - i)(A 4 - A 3 - i) 2 (A 2 - Ai - i)(\ 4 - X 1 - i)(A 4 - A 3 - i) 
11 11 



2 (A 2 - Ai - i)(A 3 - A 2 - i)(A 4 - Ai - i) 2 (A 2 - Ai - i)(A 3 - A 2 - i)(A 4 - A 3 - i) 
1 1 



(A 2 


-Ai)(A 3 


-AO(A 4 
1 


-Ai-i) 


(A 3 


-A 2 


-0(A 4 


- Ai 
1 


-i)(A 4 


- A 2 - i) 


(A 3 


-A0(A 3 


-A 2 )(A 4 
1 


- A 3 - i) 


+ (A 2 


-Ax 


-0(A 4 


-A 2 
1 


-i)(A 4 


- A 3 - i) 


(A 2 


-A0(A 3 


-A 2 )(A 4 


- A 2 - i) 


+ (A 2 


-Ai 


-<)(As 


-A 2 


-i)(A4 


- A 2 - i) 



+ 

Let us enumerate all ten terms as they enter here and make some appropriate transformations for 
each of them 



2(A 2 - Ai + i)(A 3 - Ai - i)(A 4 - A 3 + i) 

(6.11) 

1 (Ai + Q 4 (A 2 + 2)A^(A 4 + i)At 

2 (A 2 - Ai - i)(A 3 - Ax - i)(A 4 - A 3 + i) 

(6.12) 



/?> = - 




i(Ai- 


+ *) 4 (A 2 H 


-i)(A 3 + i)A|Ai 




2(A 3 


-A 2 


-i)(A4- 


Ax - i)(A 4 - A 3 


-0 


r(4) . . 
i 2 — 


i 


(Ai 


+ *) 4 (A 2 - 


hi)(A 3 + i)A|A| 




2(A 2 


-Ax 


-i)(A 4 - 


Ax - i)(A 4 - A 3 


-0 


j(4) _ 


i 


(V 


4-*) 4 (A 2 4 


-i)(\ 3 + i)\l\l 




2(A 2 


-Ax 


-0(As- 


A 2 - i)(A 4 - Ax 


-0 


J 4 — 


i 


(Ai 


+ *) 4 (A 2 - 


hi)(\ 3 + i)\ 3 3 \l 




2(A 2 


-Ax 


-0(As- 


A 2 - 0(A 4 - A 3 


-0 



2 — 2 
3—1 
4 — 3 



% (Ax + 2) 4 A3(A 3 + 2)(A 4 + «)A 4 3 



V /\ ! T ' / /x - """"" 7 ' V /x I """" ' " 
r\_/ 

i_i 2 (A 2 - Ax - i)(A 3 - Ai - i)(A 4 - A 3 - i) 

2 — 3 

3 — 4 
4—2 

(6.13) 

i (Ax + 2)(A 2 + 4 (A 3 + 2)AiAi 



(6.14) 

6 (A 2 -Ax)(A 3 -Ax)(A 4 -Ai-^) lG - i5) 



r(4 ) . (Ax + Q 4 (A 2 + z)(A 3 + 0A|Ai 
J s — i ; — rr\ ; — rr: ; 



r(4) 



(Ai+2) 4 (A 2 + ^(A 3 + 0A|A3 



(Ax + i)(A 2 + i)AlA|(A 4 + * 



(A 3 - A 2 - i)(A 4 - Ax - i)(A 4 - A 2 - i) i- 4 (A 2 - Ax - «)(A 3 - Ax - i)(A 4 - A 3 + i) 



■ 4 
2 — 1 



(6.16) 
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,(4) . (A 1 + i) 4 (A 2 + i)(A 3 + i)A|At . (A 1 + 0Af(A 2 + z)(A 3 + 4 A| 

7 l (A 3 -A 1 )(A 3 -A 2 )(A4-A 3 -i) ^ 3 '(A 2 -A 1 )(A 3 -A 1 )(A 4 -A 1 -i) [ ' ' 



r(4) i(Ai + *) 4 (A 2 + i)(A 3 + i)A|Ai i(Ai + i)(A 2 + i) 4 A|(A 4 + *)A 3 



(A 2 - Ai - i)(A 4 - A 2 - i)(A 4 - A 3 - i) ~ (A 2 - Ai + i)(\ 3 - Ai - i)(A 4 - A 3 + i) 

(6.18) 



1 2 

3^4 



.(4) (A 1+ ^) 4 (A 2 + Q(A 3 + z)A|A| (A 1+? )(A 2 + OAj(A 3 + *) 4 A| 

9 " (A 2 -A 1 )(A 3 -A 2 )(A 4 -A 2 -i) (A 2 -A 1 )(A s -A 1 )(A 4 -A 1 -0 K ' ) 



2 1 

3^2 
4^4 



r( 4) i(Ai + 4 (A 2 + i)(A 3 + i)A|Ai i(Ai + i)(A 2 + *) 4 (A 3 + i)A|A| 



(A 2 -Ai-i)(A 3 -A 2 -i)(A 4 -A 2 -i) ^ 2 (A 2 - Ai + i)(A 3 - Ai - i)(A 4 - A x - i) 

(6.20) 

Let us note that after these transformations if^ have two kinds of the denominators, namely, 
l\ A \ Jg 4 -*, 1$ have denominators of the form 

1 



(A 2 — Ai — ia)(A 3 — Ai — ife)(A 4 — Ai — id) 
with a set of integers a, b, c, while the denominators of the rest of them are of the form 

1 



(A 2 - Ai - m')(A 3 - Ai - ib')(\ 4 - A 3 - id) 

with some other set of integers a', b', d. Moreover, some of the terms 1^ have just coinciding 

denominators like, for example, and Nevertheless sometimes it will be more convenient 
to treat them separately. 

Let us start with the first group which is easier to treat. Since, the denominator of 1^ has 
the 2 <-> 3 symmetry we can simplify it as follows 



(4) . (X 1 +i)\X 2 + i)(X 3 + i){Xl + A 3 Ax + A^)A 3 
5 ~ l (A 2 -A 1 )(A 4 -A 1 -i) 



(6.21) 



where we have used the trivial identity 

A 3 = (A 3 -Ai)(A= + A3A1 + \\) + \\ (6.22) 

and the fact that the second term in the r.h.s. of (6.22) does not contribute into 
Then summing up (6.17) and (6.19) we get 

r (4) r (4) .(Ai + ^)(A 2 + 0(A2 + A : A 2 + A2)(A 3 + 4 Ai 

1-7 + in ~ 1 ~ 

7 9 (A 3 -A 1 )(A 4 -A 1 -i) 
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(Al+ , )(A2 + i) 4 (A3+ , )(A 2 + AiAg + A 2 )A 3 

(A 2 - Ai)(A 4 - 
Now adding the r.h.s. of (6.21) to (6.23) we get 

r(4) , r(4) .(4) 

. (Ai + i)(A 2 + i)(A 3 + i)((Ai + *) 2 + (Ai + i)(A 2 + i) + (A 2 + *) 2 )(A 2 + A 3 Ai + A 2 )A| 

~ -2 ; ~ 

A4 — Ai — 1 

(Ax + i){\ 2 + i)(A 3 + i)(Ai + A 2 + + A 3 Ax + \*)\\ 

~ 3 r : : = 

A4 — \\ — l 

= -3 (A 2 + i)(A 3 + i)(Ai + A 2 + *)(A 2 + A 3 Ai + + 
, q (A 2 + s)(A 3 + i)(Ai + A 2 + i)(Aj + A 3 Ax + Aj)A| 

+ 3 ; 

A4 — Ai — 1 

q (A 2 + i)(A 3 + i)(Ai + A 2 + i)(A§ + A 3 Ax + \\)\\ 

~ 6 ; ~ 

A4 — Ai — % 

q (Ai + 0(A 2 + Q(A 3 + Q 2 A 3 A 4 (A 2 + 2)(A 3 + 2)A 3 (A 2 A 3 + 1)A1 

~ 3 ; h 3 ; ~ 

A4 — Ai — 1 A4 — Ai — 1 

,1 X 2 X3 ^ o (A 2 + Q(A 3 + Q 2 A 3 A^ (A 2 + Q(A3 + 0A 3 (A 2 A3 + l)At 

~ -3A 2 A 3 A 4 + 3 ; + 6 ; 

A4 — Ai — 1 A4 — \\ — 1 



qx X 2 X s M o (A 2 + ^)(A 3 + Q 2 A 3 (A 4 + ^) 5 (A 2 + ^)(A 3 + ^)A 3 (A 2 A 3 + l)(A4 + ^) 4 

~ -dA 2 A 3 A 4 + d + 6 

A4 — Ai A4 — Ai 



1^2 

2^3 
3-^4 
4^1 



i2 x3 o ( A i + *) 5 ( A 3 + 0(A 4 + 2 A 4 (Ai + *) 4 (A 3 + i)(\ 4 + *)A 4 (A 3 A 4 + 1) 



~ — 3A 2 A 3 A 4 — 3 — 3 



A 2 — Ai A 2 — Ai 

17, x2 , 3 (± + |)(A 3 + ^)(A 4 + 2 A 4 
~ --A 2 A 3 A 4 + — + 



+ 



(3A 2 + 3iAi - ^)(A 3 + i)(A 4 + i)A 4 (A 3 A 4 + l) 17, , 2 , 3 

+ A^ To W4 + 

[t + i)(A 3 Al + 2A3 + 2a 3 A 2 -2A 2 - ? A 4 ) + (3A 2 + 3iAi - ^)(A|A| + iA 3 A| + i\\ - A 4 ) 



A 2 — Ai 

(6.24) 

In fact, it is the "canonical" form. To get this expression we have used symmetries, the formula 
(2.20) for m = 4 and m = 5 of the proposition 1 and relations (6.6) and (6.7). 
Now let us treat ij^ 

r( 4) (Ai + 0(A 2 + 2) 4 (A 3 + QAjAi . (A 1 + 2)(A 2 + Q 4 AjAj 



... 



(A 2 - Ai + i)(A 3 - Ai - i)(A 4 - Ai - i) (A 2 - Ai + i)(A 4 - Ai - i) 
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.(\ 2 + i) 4 \ 3 3 \l , (A 2 + *) 4 A 3 A 4 £ (A 2 + *) 4 A 3 A 4 



A 2 - Ai + i (A 2 - Ai + i)(A 4 - Ai - i) (A 2 - \\ + i)(A 4 - A x - i) 

, A 4 A 3 (A 4 + Q 4 , A^(A 4 A 4 + 4a 3 A 4 - 6A^A 4 - 4?A 2 A 4 + A 4 ) 

'(A 2 -A 1 )(A 4 -A 1 ) ~ 1 (A 2 -A 1 )(A 4 -A 1 ) 

.A 3 (4*A 3 A 3 - + Ax) - 4^A 2 A 4 (A2 + A 4 Ax + A?) + A 3 + A4A1 + A 4 A 2 + A?) 

A 2 — Ai 

.Ai(-6AiA|A| - 4iAiA 2 A^ - 4iA 2 A 2 A 4 + AiAf + A 2 A 4 + A?) 

A 2 — Ai 

iA|j(-6A?A;{ - 4iAiA^ - 4iA 2 A 4 ) + i*'' 



2x2 ^ x2 . . • Ai(-6AfA| - Ai\\\\ - Ai\\\ A + A X A| + A 2 A 4 + A?) 



A 2 — Ai 

^ ^ A 3 (A 2 (5^A 2 - 5Ai - %i) + A 4 (-5A 2 + 4zAi + 1) - fiAj + Ai + |) 

A 2 — Ax 



+ A^(-3)A^A 2 + iA 3 (-4)iA 4 ^A 2 



3 (5A 2 + gAi - |)A|A| + (5*A 2 - 4Ai - *)A 3 A 3 + (-|Aj - ^ + |)A| 

~ — A 2 A 3 A 4 + [O.ZO) 

A 2 — Ai 

Here the symmetry, the formula (2.20) for m = 4 and the relation (6.1) were used. 

Now we shall treat the other six terms /j 4 -*, J 4 4 ' ) , 1^ and 1^ given by the expressions 

(6. 11), (6. 12), (6. 13), (6. 14), (6.16) and (6.18) respectively. Here we shall proceed in two steps. As 
a first step we will reduce all these six terms to a following form 

r(4) A? Cj Dj 

3 ~ (A 2 -A 1 )(A 3 -A 1 )(A 4 -A 3 ) + (A 2 -A 1 )(A 3 -A 1 ) + (A 2 - Ai)(A 4 - A 3 ) + A 2 - \ 

(6.26) 

where j = 1, 2, 3, 4, 6, 8 and Aj, Bj, Cj, Dj are some polynomials. Then we shall sum up all the six 
results and get the "canonical" form for the sum obtained. 
Let us start with the expression (6.11) for the term 

r( 4) i (Ai + i)A?(A 2 + i)Ai(A 4 + i) 4 



2 (A 2 - Ai + i)(As - Ai - i)(A 4 - A 3 + i) 

j (Aj + Ai(A 2 + z) + (A 2 + if)(\ 2 + QA|(A 4 + Q 4 
2 ' A, - A ;! • / 



+ 



1 (A 2 + Ai(A 2 + i) + (A 2 + i) 2 )(A 2 + i)A 4 (A 4 + i) 4 i (A 2 + i) 4 A 3 (A 4 + i) 3 

2 (A 3 - Ai -i)(A 4 - A 3 + i) 2 A 2 - Ai + % 
i (A 2 + i) 4 A 4 (A 4 + i) 3 i (Ai +i)(A 2 + i) 4 A 4 (A 4 + i) 3 



2 (A 2 - Ai - i)(A 3 - Ai + i) 2 (A 2 - Ai + i)(A 3 - Ai - i)(A 4 - A 3 + i) 



i (\\ + Ai(A 2 + i) + (A 2 + i) 2 )(A 2 + i)A|A 



3\4 
4 



2 A 4 — A 3 

1 (Af + Ai(A 2 + 1) + (A 2 + Q 2 )(A 2 + i)(\ 3 + Q 4 (A 4 + i) A i A 4 A|(A 4 + ?) 3 

2 (A 3 -Ai)(A4-A 3 ) 2 A2-A1 
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+ 



_ i A 4 (A 3 + ?) 4 (A 4 + ?) 3 i {\i+i)\\{\z + i)\K + if 
2 (A 2 - A0(A3 - Ai) 2 (A 2 - Ax)(A 3 - Ax)(A 4 - A 3 ) ~ 

1 A 3 A 4 (A 2 + A 3 (A 4 + i) + (A 4 + *) 2 )(A 4 + i) 

r^j -j- 

2 A 2 — Ai 

j (Ai + ?) 4 (A 2 + Q 4 (A| + A 3 (A 4 + 1) + (A 4 + Q 2 )(A 4 + i) _ i A 4 Aj(A 4 + z) 3 
2 (A 2 -A 1 )(A 3 -A 1 ) 2 A 2 -Ai 

^ A 4 (A 3 + *) 4 (A 4 + *) 3 ^ (Ai + *)A 4 (A 3 + *) 4 (A 4 + tf 



2 (A 2 - A0(A 3 - Ai) 2 (A 2 - Ai)(A 3 - Ai)(A 4 - A 3 ) 

where we have used symmetries and the formula (2.20) with m — 4. 
For the next term we start with the expression (6.12) 



7 (4) j (A 1 +i) 4 (A 2 + i)Aj(A 4 + i)Ai 
2 ~ 2(A 2 -A 1 -i)(A 3 -A 1 -i)(A 4 -A 3 + i) 

i (Ai + 4 (A 2 + *)(A 2 + A 3 (A 4 + i) + (A 4 + *) 2 )(A 4 + i)A^ 



+ 



i (Ax + z) 4 A 3 (A 4 + z)A 3 i (A : + z) 4 A 3 (A 4 + z)A 4 



4 



(6.27) 



2 (A 2 -Ai-i)(A 3 -Ai-i) 
i (Ax + 2) 4 (A 2 + 2)(A 4 + 2) 4 A 3 
+ 2 (A 2 - Ai - i)(A 3 - Ai - z)(A 4 - A 3 + i) 

i A 4 (A 2 + Q(A 2 + A 3 (A 4 + i) + (A 4 + Q 2 )(A 4 + i)A| % A 4 (A 2 + QA 4 (A 4 - if 
~2 (A 2 -A 1 )(A 3 -A 1 ) 2(A 2 -A 1 )(A 3 -A 1 )(A 4 -A 3 ) l ' j 

The term given by (6.13) is also simple to treat 

7 ( 4 ) = j (Ax + i) 4 A|(A 3 + i)(A 4 + i)A| 
3 2 (A 2 - Ai - i)(A 3 - A x - i)(A4 - A 3 - i) 



2 (A 2 - Ai - i)(A 3 - Ai - i) 2 (A 2 - Ax - i)(A 3 - Ax - z)(A 4 - A 3 - i) 

z A 4 Aj(A 4 + QA| i AfA|(A 4 + 2z)(A 4 + z) 4 
~ 2(A 2 -A 1 )(A 3 -A 1 ) + 2(A 2 -A 1 )(A 3 -A 1 )(A 4 -A 3 ) 1 ' J 

The fourth term tf^ (6.14) demands more work in order to reduce it to the form (6.26) 



J 4 ~ 



1 (A 1+ 2)(A 2 + z) 4 (A 3 + z)A|A 3 

2 (A 2 - Ax + i)(A 3 - Ax - i)(A 4 - A 3 - i) 



(Ax+2)(A 2 + 4 A|A 3 i (Ax + *)(A 2 + *) 4 A 3 A 4 



+ 



2 (A 2 - Ax + i)(\ 3 - Ax - 2) 2 (A 2 - Ax + i)(A 3 - Ax - i)(A 4 - A 3 - i) 
^(A 2 + .) 4 A 3 A 3 i (A 2 + .) 4 A 4 A| 



2(A 2 -Ax + z) 2(A 2 -Ai + z)(A 3 -Ai-z) 



+ 



i(Ai + i)(A 2 + *) 4 (A 2 + A 3 (Ax + i) + (Ax + *) 2 )A 4 , i(Ai + *) 4 (A 2 + *) 4 A 4 



+ 



2(A 2 - Ax + i)(A 4 - A 3 - i) 2(A 2 - Ax + i)(A 3 - Ax - i)(A 4 - A 3 - i) 
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i A 4 (A 3 + if\\ i (Ai + i)A|(A| + A 3 (Ai + i) + (Ai + z) 2 )(A 4 + z) 4 

~ -7:7^ t^t^ r^r + t:- 



2(A 2 -A 1 )(A 3 -A 1 ) 2 (A 2 -A 1 )(A 4 -A 3 ) 

_ z A 4 A 2 (A 4 + z) 4 

~ 2(A 2 -A 1 + 0(A 3 -A 1 )(A 4 -A 3 ) ~ 
i Aj(A 3 + i) 4 Aj z (Ai + z)Af(A| + A 3 (Ax + + (A x + z) 2 )(A 4 + if 

2(A 2 -A 1 )(A 3 -A 1 ) + 2 (A 2 -Ai)(A 4 -A 3 ) 
i(A? + A 2 (A 2 + z) + Ai(A 2 + z) 2 + (A 2 + z) 3 )A 4 (A 4 + z) 4 z(A 2 + ifX 4 2 (X 4 + tf 



2(A 3 - AO(A 4 - A 3 ) 2(A 2 - Ai + z')(A 3 - Ai)(A 4 - A3) 

^ A 4 (A 3 + z) 4 A 3 z (A 1 +^)A 4 (A 2 + A 3 (A 1 +0 + (A 1 +2) 2 )(A 4 + 4 



2(A 2 -A 1 )(A 3 -A 1 ) 2 (A 2 -AO(A 4 -A 3 ) 
i (A 2 + z) 4 A 4 (A 3 + A 2 (A 4 + z) + A3(A 4 + z) 2 + (A 4 + z) 3 ) z A 4 (A 2 - z) 4 (A 4 + if 



2 (A 2 -A 1 )(A 3 -A 1 ) 2(A 2 -A 1 )(A 3 -A 1 )(A 4 -A 3 ) 

(6.30) 

Now we treat ijf^ given by (6.16) as follows 

(4) (\ 1 +l)(\ 2 +l)\ 3 2 \l(\ i + lf 

6 ~ *(A 2 -A 1 -z)(A 3 -A 1 -z)(A 4 -A 3 + z) 

= . (Aj + z)(A 2 + z)(A 2 + A 2 (Ai + z) + (Ai + z) 2 )A|(A 4 + if 
(A 3 - Ai -z)(A 4 - A 3 + z) 

(A!+2) 4 (A 2 +0A3(A 4 + 2) 4 



(A 2 - Ai - z)(A 3 - Ai - z)(A 4 - A 3 + z) 

, (A 2 + z)(A 2 + A 2 (Ai + z) + (Ai + z) 2 )A 3 (A 4 + z) 4 
A 4 - A 3 + z 

, . (A 2 + z)(A 2 + A 2 (A! + z) + (Ax + z) 2 )A 4 (A 4 + if (Ai + z) 4 (A 2 + z)A 3 (A 4 + z)' 



(A 3 - Ai - z)(A 4 - A 3 + z) (A 2 - Ai - z)(A 3 - Ai - z)(A 4 - A 3 + z) 

= , (A 2 + z)(A 2 + A 2 (Ax + z) + (Ax + z) 2 )A|A 4 | 
A 4 — A 3 

, . (A 2 + z)(A 2 + A 2 (Ai + z) + (Ai + z) 2 )(A 3 + z) 4 (A 4 + z) 4 . A 4 (A 2 + z)A|A 4 

\ 1 /\\\/\\\ \ 1~ 



(A 3 - Ax)(A 4 - A 3 ) (A 2 - Ax) (A 3 - Ax)(A 4 - A 3 ) 

_ , A| A 4 (A 3 + z)(A| + A 3 (A 4 + z) + (A 4 + z) 2 ) 

A 2 — Ai 

_ . (Ai + z) 4 (A 2 + z) 4 (A 2 + A 4 (A 3 + z) + (A 3 + z) 2 )(A 4 + z) A 4 (A 2 + z)A 3 A 4 

(A 2 -A 1 )(A 3 -A 1 ) +\a 2 -A 1 )(A 3 -A 1 )(A 4 -A 3 ) 1 • j 

The last term 1^ (6.18) is more simple 

r( 4) . (A 1 +z)(A 2 + z) 4 A 3 (A 4 + z)Ai 

Jo ~ —%- 



(A 2 - Ai + z')(A 3 - Ai - z)(A 4 - A 3 + z) 
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. (A 2 + i) 4 A|(A 4 + i)Ai . (A 2 + *) 4 A 4 (A 4 + *)Ai 



(A 2 - Ai + i)(A 4 - A 3 + i) (A 2 - A x + i)(X 3 - X 1 - i)(X 4 - A 3 + i) 
. (\ 2 + iy\ 3 3 \l (A 2 + i) 4 AjA| (A 2 + *) 4 A!(A 4 + *)Ai 



A 2 -Ai+i (A 2 - Ai + i)(X 4 - A 3 + i) (A 2 - Ai + i){X 3 - Ai - i)(A 4 - A 3 + i) 

■ A 4 (A 3 + if X\ . A 4 (A 3 + *) 4 (A 4 + *)A| 

~ *(A 2 -Ai)(A4-As) '(A 2 -A 1 )(A 3 -A 1 )(A 4 -A 3 ) ^ d j 

Now we are prepared to perform the next our step. Namely, we will gather all the six results 
(6.27-6.32) into the form like (6.26) 

4 4) + 4 4) + 4 4) + 4 4) + 4 4) + 4 4) ~ 

A B C D 

+ 7^ \~TT\ TT + 7^ \~TT\ TT + 1 r (6-33) 



(A 2 -A 1 )(A 3 -A 1 )(A 4 -A 3 ) (A 2 -A 1 )(A 3 -A 1 ) (A 2 - A X )(A 4 - A 3 ) A 2 - X l 
where 

A = - i -(X l +i)Xt(X 3 + t)\X A + if + 1a 4 (A 2 + *)A 4 (A 4 -*) 3 + i -XtX 3 (X 4 + 2t)(X 4 + tf + 

+ ^ 4 (A 2 - i)(A 4 - if + «A 4 (A 2 + *)A 3 A 4 - *A 4 (A 3 + *) 4 (A 4 + i)X\ (6.34) 

B = ^(Ai + i)\X 2 + *) 4 (A 2 + A 3 (A 4 + i) + (A 4 + *) 2 )(A 4 + i) - *A 4 (A 3 + *) 4 (A 4 + if + 

+ ^ 4 (A 2 + i)(A 2 + A 3 (A 4 + i) + (A 4 + *) 2 )(A 4 + *)A 3 + *A 4 A 3 (A 4 + i)X\ + 

+ i(Ai + *) 4 (A 2 + *) 4 (A 2 + A 4 (A 3 + i) + (A 3 + *) 2 )(A 4 + i) (6.35) 
C = ^(Ai + t)Xt(Xl + A 3 (Ai + i) + (Ai + *) 2 )(A 4 + if + *A 4 (A 3 + ifX\ (6.36) 

D = t -X 3 1 Xj(X 2 3 + X 3 (X 4 + i) + (X 4 + if)(X 4 + i) - *A 4 A 3 (A 4 + *) 3 - 

- *A 3 A 4 (A 2 + A 3 (A 4 + 1) + (A 4 + *) 2 )(A 3 + i) (6.37) 

Now we want to get the "canonical" form (2.7) for the expression (6.33). To do this we actively- 
used the program MATHEMATICA because the calculations are straightforward but become more 
cumbersome. Let us outline our further actions. 

It is more convenient to start with the first term in the r.h.s. of the formula (6.33) 

A 

(A 2 -A 1 )(A 3 -A 1 )(A 4 -A 3 ) (6 ' 38) 
We shall use the fact that the denominator is antisymmetric under the substitution 

Ai <-> A 3 
A 2 <-> A 4 
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Since A is a polynomial given by (6.34) then the simplification procedure of the term (6.38) is as 
follows. If in the expression (6.34) one faces a monomial 

A^A^AfAl 4 (6.39) 

where without loss of generality 24 > i 2 one can apply the evident identity 

{t\ \ \ v^«4-i2-l \ i'2+k \ i 4 -i 2 -l-fe • v. „• . 

{A4 — A 3 ) 2^k=0 A 4 A 3 1 M > *2j (6.40) 

, *4 = *2 

Therefore if i 4 > i 2 then 

\«1 \«2 \«3 \M \il \12 \M— 12+13 \ 12 *2~ 1 \ U \ *2 \ «4 — *2+«3 — 1 — \ *2+& 

Aj A 2 A ;j A j A , A 2 A ;j Aj ^ Aj AoA :j A , 



(A 2 - Ai)(A 3 - Ax)(A 4 - A 3 ) (A 2 - Ai)(A 3 - A X )(A 4 - A 3 ) fo ~ A i)( A 3 - Ai) 

(6.41) 

Let us note that the second term in (6.41) gives rise into the second term " B" in the r.h.s. of the 
formula (6.33). If i 4 = % 2 then only the first term in (6.41) survives. The first term in (6.41) is 
symmetric under the transposition A 2 <-> A 4 . If i± + i 2 = H + ii then it is also symmetric under the 
transposition Ai <-> A3 and the first term is "weakly" equivalent to zero according to the formula 
(2.15). If i 1 + i 2 < «3 + 14 the following "weak" equality is valid 



\ U \ 12 \14— «2+«3 \*2 1 / \ n \ «3+*4-«2 \W \ I3+I4 — »2\ \«2 \« 

A^ A 2 A3 A4 1 ^A^ A3 — A3 A^ / / *2 A4 



(A 2 - A X )(A 3 - A0(A 4 - A 3 ) 2 (A 2 - Ai)(A 3 - A X )(A 4 - A 3 ) 

Y l3+*4— 12— »1— 1 ^«l+fc^i2 ^«3+«4— «2 — 1 — fc^l2 

= 2 £ 1 (A 2 -A 1 )(A 4 -A 3 ) 4 (6 ' 42) 
For the case i\ + i 2 > i 3 + 14 the sum in (6.42) should be substituted by 

1 ii+i 2 — «3~ *4 — 1 \U— 1— fc \»2 \«3+24— 12+k \ 12 

-- T ^ ^ ± (6.43) 

2 to (A 2 -A 1 )(A 4 -A 3 ) 1 ] 

In both cases, namely, if i\ + i 2 ^ iz + 14 the sum in (6.42) or (6.43) gives rise into the third term 
"C" in the r.h.s. of (6.33). 

Performing this procedure for all monomials of the form (6.39) participating in the polynomial 
A given by the formula (6.34) one can arrive at the formula 

A (4) + 4 4) + 4 4) + if + 4 4) + 4 4) - 7, tt + 77 rr + D 



(A 2 - Ai)(A 3 - AO (A 2 - A0(A 4 - A 3 ) A 2 - X 1 

(6.44) 

with some other polynomials B' and C . 

Due to the 2^3 symmetry of the denominator of the first term in (6.44) one can treat any 
monomial A^ 1 A2 2 A3 3 A 4 4 participating in B' as follows 

\jl\j2\j3\j4 33-32-1 \ 31 \ 32+k \ 33-k-l X 3\ 

A 1 A 2 A 3 A 4 ^ A 1 A 2 A 3 A 4 . . 

L x x ( 6 - 45 ) 



(A 2 — Ai)(A 3 — Ai) k=0 A 2 — Ai 

where without loss of generality it is implied that j 3 > j 2 because if j 2 = j 3 then due to the 2^3 
symmetry and the formula (2.15) this term would make zero contribution. So the r.h.s. of the 
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(6.45) gives rise into the third term in (6.44). Using this one can treat the whole first term in 
(6.44). 

For any monomial A^A^A^A^ 4 participating in C of the second term of the expression (6.44) 
one can write 

\k1\k2\k3\k4 k-2 — 1 \ ki+l \ k2 \ k:j — 1— I \ £14 k±— 1 \ fci+fo \ ks+l \ k^ — 1— I \ ki+k2 \ £3+^4 

A l A 2 A 3 A 4 A l A 2 A 3 A 4 . \p A l A 3 A 4 A l A 3 



(A2 — Ai)(A4 — A3) i=Q A4 — A3 i=Q A2 — Ai (A2 — Ai)(A4 — A3) 

(6.46) 

fe2— 1 \ fcl \ k4 \ ki+l \ k2 — 1— / k^ — 1 \ fci+fc2 \ fe3+/ \ ^4 — 1— i \ fcl+/C2 \ &3+fc4 
„ £ Al A2 A 3 A 4 + £ ^1 A 3 A 4 + Ai A3 (g 4 ^ 

; =0 A2 — Ai ;=0 A2 — Ai (A2 — Ai)(A 4 — A3) 

where for the first sum of (6.46) we have applied transformation 

Ai <-> A 3 
A 2 <-> A 4 

So the both the first term and the second term in (6.47) give rise into the third "D" term in (6.44) 
while the third term in (6.47) gives rise to the second term "C" in (6.44). Proceeding this way 
one can treat the whole expression 

a 



(A 2 -Ai)(A 4 -A 3 )' 
As a result of performing this scheme one can arrive at the expression 

if + if + if + if + 4 4) + if - 

C" D' 



(6.48) 



(A2 — Ai)(A4 — A3) A2 — Ai 

where C" is a polynomial of two variable Ai and A3. 
Now with the help of the identity 

— \h+i2 _|_ I (A2 - Ai) EfcLo 1 Ai 1+fc A 2 2 1 k , %2 > 0; (6.49) 
one can reduce the second term in (6.48) to the form 

jy d" 



A2 — Ai A2 — Ai 



+ E 



where D" is a polynomial of Ai, A3 and A 4 and E is some polynomial. 

What is left now is to reduce the power of the polynomials C", D" and E with the help of 
the formulae (2.20), (2.21), (2.22) or (6.1-6.4). We should also use the fact that if we do the 
substitution Aj — > Xj — i/2 there is a restriction that the function should be even under the 
transformation {x±, x 2 , x 3 , X4} — > {— x±, — x 2 , — x 3 , — X4} according to the formulae (2.9), (2.10) 
and (2.11). 

As a result of all these actions described above we get the " canonical" form for the sum 

if + If + If + If + If + If ~ C '" - + ^ + Bf (6.50) 

{A2 — M){M — A3J A2 — Ai 
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where 



where 

C" = 2XlX 2 3 + AtX 1 X 2 3 - ^X 2 3 - ^X 4 X 3 - iX 3 + 1 (6.51) 

29 Fi 
D'" = Xl(22XjXl + 22iX 3 X 3 4 - yA 3 + 19X 3 X 2 4 - °-iX\) + 

, ^.,2,^ 47, ,s 61. ,o 67., , 2 3, 2 23., , 
+ A x {22iX\X\ - — A 3 A 3 - — *A 3 + — iX 3 X\ - -X\ + — tX 4 ) - 

88, 2 , 3 367., ,o 427, o 303 , , 2 37. , 2 97, , ,„ rrt , 

AoA? tX 3 X 3 4 + A4 A3A4 + — iXl A 4 6.52 

5 3 4 20 4 40 4 40 4 40 4 20 ; V ; 

E> = - — A2A3A4 (6.53) 

Now we have to sum up the four contributions (6.10), (6.24), (6.25) and (6.50) and get the 
"canonical" form for T 4 (Ai, A 2 , A 3 , A 4 ) 

p(4) p(4) 

T 4 (A 1 ,A 2 ,A 3 ,A 4 ) ~ P (4) + + n TTa \~T (6-54) 

P (4) = -f A^A 3 (6.55) 

?| 4) = + 30 iX 3 X 3 4 - 16XI + 18 A^ + 8A 4 ) + 

+ Ai(30iA^A 3 + 30A 3 A 3 - 16 iX\ + 18iA 3 A^ - 4A^ + 4iA 4 ) - 

^4 42 43 

-20A^A 3 - 20*A 3 A 3 + —A 3 - — A 3 A^ - — *A 4 (6.56) 

P 2 {4) = 2X\X\ + AiX^l - ^X 2 3 - ^X 4 X 3 - tX 3 + i (6.57) 

Let us note that in terms of the real variables Xj the polynomials (see (2.11)) look a little 
bit simpler 

o 

P{ = x 1 (30^3X4 - — x'l + —x 3 x 4 + —x 4 ) - 4xix 4 - 

2 2 8 

25 2 o 147 . 531 o 653 
- y« + ~^x 4 - ~^x 3 x A - — x 4 (6.59) 

P 2 (4) = 2a^ + X - Xl x 3 -lx 2 3 + ^ (6-60) 



So, the function 

>(4) p(4) 
+ — 

x 2 -x 1 (x 2 - xi)(x 4 - x 3 ) 

is even i.e. 



P l 

T 4 (xi,x 2 ,a; 3 ,a; 4 ) = P (4) + + £ (6.61) 



T 4 c (-xi, -x 2 , -x 3 , -x 4 ) = T 4 c (x 1 ,x 2 , £3,2:4) 
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as it should be according to the formula (2.13) with n — 4. 

Now let us start the second step of our general scheme, namely, the performing of the integra- 
tion of the "canonical" form (6.54) 



p ( 4 ) = II / ^^4(A 1 ,A 2 ,A 3 ,A4)T 4 (A 1 ,A 2 ,A3,A 4 ) = J (4) + j{ 4) + ^ 4) (6.62) 

j = l J C-l/2 



where 

4 d\ 



4 4) = II / P: ^(Ai, A 2 , A 3 , A 4 ) P (4) (6.63) 

4 p(4) 

4 4) = II / ^ ^(Ai, A 2 , A 3 , A 4 ) (6.64) 

4 J\ p(4) 

J 4) = II / ^(Ai, A 2 , As, A 4 ) — (6.65) 

fJ[Jc_ 1/2 2m (A 2 - Ai)(A 4 - A 3 ) 

Using the formulae (6.63-6.65), (6.55-6.57), (2.46), (2.47), (2.55) we get 

4 4) = II / ^ ^(Ax, A 2 , A 3 , A 4 ) (-^)A 2 A^ = 
jJ[Jc_ 1/2 2m 5 

OO OO OO OO OA 

= £ (4) E ("I)" E ("I)' 2 E ("I)' 3 E (-l) U (-^)(^2 + e 2 )(ih + e 3 ) 2 (^4 + e 4 ) 2 = 

Z 1= Z 2 =0 i 3 =0 k=0 

oo oo oo oo q 4 

= ^EC-i^El-i^EC-i^EC-i)^-^ 2 ^ = 

; 1= o i 2 =o h=o h=o 

1 -p-r / d d x . 34, 2 3 17 ,„ „„, 

Restoring for convenience the dependence of the polynomials P^ and on A-s as in the 
formula (2.8) we can get for the term 

d\ rT „ , x , , A (4) (A 1 |A 3 ,A 4 ) 



= n I Pl Ua{XuX2jX37Xa) 

i J C_ i /o 



3 JiJc_ 1/2 2m A 2 - Ai 

= j>2( i) h E( 1 )' 2 E( 1 )' 3 E( l)' 4 - ^ 4 ^ - ^ 1 + - tls + 63 '^ 4 + e ^ = 

h=i i 2 =o i 3 =i « 4 =o ^2 + iii + e 2 - ei 

/•l /^Q oo oo oo oo 

6 h=l «2=0 J 3 =l k=0 

P^i-ih + e 1 \-il 3 + e 3 ,tl 4 + e 4 )s' 1+i2+l(ei - £2) = 
/■i/ 15(s-l) 3 (7 - 26s + 7s 2 ) , 21 (s - 1) , 2 7 (2 - s + 2s 2 ) , 3 
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+ 5 lnS + 1 (1 + 223+^ > = 5 61 _ 65 

48(l + s) 3 240 (1 + s) 4 > 8 20 s v ; 32 s v ; v ; 

The last term ji 4; can be calculated as follows 



4 4) = n / ^ ^(A 1? A 2 , As, A 4 ) - P * (Al ' A 

1 J C_i /o Z7TZ 



3 J 



fJ[Jc_ 1/2 2ni ' ' ' (A 2 - Ai)(A 4 - A 3 ) 



>(*)/ 



= ^ (4) E(-i) h E(-i)' 2 E(-i)' 3 E(-i) S 7 t fv 1 \r~l 3 r 1 3) \ = 

i% h^o h^i i% {il2 + ih + e2-e 1 )(il 4 + il 3 + e 4 -e 3 ) 

s - 70 1 h=l l 2 =0 h=l l 4 =0 

P 2 (4) (-iZi + ei, -iZ 3 + e 3) s i i+ / 2^(^-^) t '3+/4+i( £ 3-e4) = 

= /^ds / 1 rft(4 / lnS( ' ln J ~ ln ^ + (--(-3-s + t + 3st)ln 3 s + 4(-l + 3s - 3t + st) In 2 s hit 
jo jo L (l + s) 2 (l+t) 2 3 

3 (3- 13s + 3t + 3st)ln 2 sln 2 t + -^-(3 - 5s - 5t + 3st) In 4 s)- 



4 V 7 12 v 7 7 (l + s) 3 (l + t) 3 

1 (-1 + 5s - 2s 2 + 2t- 5st + s 2 t) , „ 2 , 2 , 2 _ , 2 x 

- In s(ln 2 s - In 2 1) (In 2 s - 5 In 2 1) + 



6 (l + s) 4 (l + t) 3 

+ (-(3 + 30s - 5s 2 - 34t - Ast - UsH - 5t 2 + 30st 2 + 3s 2 t 2 ) In 3 sin t + 
3 

+ — (4 - 17s + 9s 2 - lit + 66st - 17s 2 t + 9t 2 - 17st 2 + 4sV) In 5 s hit - 

- — (2 - 31s + 7s 2 + Ut + 33st - 21s 2 t + 2t 2 + 4st 2 + 2s 2 t 2 ) In 3 s ln 3 f)- -} 

15 v ' 7 (1 + s) 4 (l + 1) 4 J 

Taking the integral we come to an answer for J 2 4; 

4 4) = 4C(3) - yC(3)ln2 - |C(3) 2 - |C(5) + |c(5)ln2 (6.68) 
Summing up all the three results (6.66), (6.67) and (6.68) we finally get our main formula (1.9) 

p(4) = jr + j; 4) + 4 4) = 

1 173 11 51 55 85 

L _ 21n2 + — C(3) - — C(3) ln2 - — ( 2 (3) - — C(5) + — C(5) In 2 (6.69) 

5 60 v ; 6 v ; 80 v ; 24 sv ; 24 sv ; v ; 
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